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Bose gas in power-like spherically symmetric potential in arbitrary spatial 

dimensionality. 
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Temperature of the Bose - Einstein condensation and the temperature behavior of the chemical 
potential and other thermodynamical functions of the ideal Bose gas are found for the arbitrary 
power-like spherical-symmetric potential at an arbitrary space dimension. It is shown that the 
recently observed Bose - Einstein condensation of photons in the cavity is the phase transition of 
the third kind. 
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I. INTRODUCTION 



The majority of studies, both theoretical and experi- 
mental, on the Bose - Einstein condensation (BEC) was 
carried out for gases in the harmonic traps and in three 
space dimensions [l|. In the meantime, a number of 
works were devoted to the study of the thermodynamic 
properties of the Bose gas trapped in arbitrary power- 
like spherical-symmetric potential in three spacial (3D) 
dimensions |2j, and to the harmonic trap at lower dimen- 
sion D = 2 Q. Recent publication [J] contains the cal- 
culation of the temperature and the fraction of the con- 
densed particles in case of Bose - Einstein condensation 
of the ideal Bose gas in free space at any space dimen- 
sion. As far as the physical realization is concerned, the 
harmonic traps of lower dimensions are effectively mod- 
eled by the sharply anisotropic frequencies 0]. A direct 
realization of BEC in the 2D harmonic and isotropic trap 
was recently achieved with the photons confined in the 
spherical cavity Q. The confinement results both in the 
effective nonrelativistic dependence of the photon energy 
on the transverse momentum, the harmonic oscillator po- 
tential in the transverse direction, 
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with a tiny effective photon mass m e g 
5 x 1CP 33 gramm, the frequency lj — C\j2j (RD ) ~ 
4 x 10 11 s" 1 , and E = nhcs/D « 3 eV |. Here, 
Dq = 1.4 x 10 ~ 4 cm is the distance along the optical 
axes between the spherical mirrors, each having the ra- 
dius of curvature R — 100 cm, s = 7 is the excitation 



number of the longitudinal standing wave of the photon, 
and r, k r = p r /h are, respectively, the distance from the 
optical axes and the conjugate wave number in the trans- 
verse direction, c and Ti are the velocity of light and the 
Planck constant. 

However, the general analysis of the possibility of BEC 
in a power-law potential in an arbitrary space dimension 
D, as well as the degree of singularity of the chemical 
potential and other thermodynamic functions, is absent 
in the literature. In the present work we fill the gap and 
find both the conditions for BEC, as well as the degree of 
singularity of the thermodynamical functions of the ideal 
Bose gas, appearing in the case of arbitrary power-like 
spherical potential and for arbitrary space dimensional- 
ity. Despite the fact that the ideal Bose gas approxima- 
tion adopted in the present work is poor for the gases of 
alkali metals, it is meaningful for BEC of photons stud- 
ied in Ref. [f| , in addition to the case of the hypothetical 
weakly interacting axions, proposed as dark-matter can- 
didate in Ref. @. 

The rest of the paper is organized as follows. The den- 
sity of states is obtained in Section |TTJ Sec. Mil is devoted 
to the elucidation of the behavior of the chemical poten- 
tial near the the BEC transition temperature. The re- 
sults and conclusions are presented in Sec. [IVJ Appendix 
contains the necessary information concerning the poly- 
logarithm function. 



II. DENSITY OF STATES AND THE CRITICAL 
TEMPERATURE OF BEC 



The starting point is the semi-classical one-particle 
density of states v(e) in the case of spherically symmetric 
single-particle potential of the form 



(2.1) 
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characterized by its strength Uq > 0, the scale a, and the 
exponent 7. In effect, the single-particle density of states 
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can be denned as 

d D rd D p , 



u{e) 



D 



P 

2m 



U, 



■•(;)' 



(2.2) 



where S(x) is the Dirac delta function. Hereafter D is 
the space dimensionality of the gas, the particle spin is 
assumed to be zero, and p and r are the moduli of the 
_D-dimensional vectors of the momentum and the radius, 
respectively. The phase-space integral (|2.2I) is calculated 
in the multi-dimensional spherical coordinates, where the 
volume element integrated over angular variables takes 
the form 
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See 0, the formula 4.632(2) differentiated over R. Anal- 
ogously for d D p. With the use of the table integral [|| 
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the calculation leads to the explicit expression 
v(e) = At?- 1 , 



(2.3) 
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and r(x) is the Euler gamma function. 

The case of 7 = 0, at first sight, seems to be singular. 
However, it corresponds to the situation, when the gas 
is under the influence of the constant potential Uq. As 
is usual in the case of free particles or particles in the 
constant potential, one should put the system into the 
finite volume. The problem is reduced to the case of 
bosons in D spatial dimensions, with the only difference 
that they condense to the state with the energy Uq. The 
temperature of BEC in this case is the same as found in 
4]. The thermodynamical properties of the system do 
not depend on changing the energy reference point. 

In the meantime, there is another way to get the case 
of free bosons in the space with no external potential, 
confined inside the sphere with the radius a. This is 
obtained in the limit 7 — > 00. Then 



U(r) = U (^ 



0, if r < a, 
00, if r > a, 



(2.5) 



The density of states in this limit coincides with that 
obtained earlier [4] for free bosons in a finite volume. 
As for the dependence of the chemical potential on the 
temperature, in the cases of 7 = or 7 — > 00, we will 
comment on this point after establishing the analogous 
dependence in the general case of finite 7 ^ 0. 

As usual, the inequalities n(T) < and (fy) < 0, 

result from the condition of having the particle number 
bosons in the gas, i.e., 
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where p, = fi(T) is the chemical potential. The above 
equation is valid only if Uq > 0. The Boltzmann constant 
is set to unity. BEC is possible if the integral defining 
the temperature of BEC Tq, 
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has a definite value. Hereafter 



00 1 



is the Riemann zeta function. This function is defined at 
x > 1. Hence, the restriction on the parameters of the 
problem 7 and D, resulting from the definition of zeta 
function, is 



D D , 



(2.9) 



One obtains that BEC is possible in the space dimension 
D = 1,2, and 3 at, respectively, 7 < 2, 7 > 0, and 
7 > —6. In particular, the confinement in the harmonic 
trap, 7 = 2, admits BEC in any space dimension when 
D > 1 . The critical temperature T of BEC is found from 
Eq. 
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In particular, one has D — 3, 7 = 2 for BEC of alkali 
atoms in the traps with the typical magnitude of the 
frequency uj — 10 2 s^ 1 . Then ma 2 /(2h 2 U a ) = {huj)~ 2 , 



and the expressions (|2.4[) are reduced, respectively, to the and the BEC temperature is 
following ones 
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for N ~ 10 6 . Hereafter, the necessary expressions for the 
zeta function at particular arguments are taken from [9[ . 
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On the other hand, under the conditions of the experi- 



2,7 = 2, 



10 1 



mcnt [5[ , the parameters are D 
and the BEC temperature is 



Tq = — V3N - 10 J K, 

at N ~ 10 4 photons in the cavity. The above equation 
allows for two polarization states of the photon. For the 
oscillator potential in space with arbitrary dimension D 
one finds 

l/D 

(2.11) 
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Since lirno-^ooCP) = 1) one can see from (|2.11|) that, in 
this limit, To — > %uj independently of the particle number 
N. 

The chemical potential n(T) at T < To is zero with all 
its derivatives. Hence, the number of particles sitting at 
the excited energy levels, JV>(T), is given by Eq. (|2.8[) . 
in which one should make the replacements N — > N> (T) 
and To T. Then, the number of particles in BEC, 
Nq(T) = N — iV>(T), is given by the expression 



iV (T) = N 



(2-12) 



This expression is valid for the finite values of the power 
7 in the external potential. It coincides with the earlier 
expression a = 3/2 + 3/7 [H,[i] at D = 3. As is explained 
above, the case of ideal Bose gas put in a finite volume, 
upon neglecting the interaction with the external field, 
corresponds to taking the limit 7 — ¥ 00, and the power a 
is reduced to the value a — 3/2 in the three-dimensional 
case. 



Excluding from these equations the multiplicative factor 
A, one obtains the internal energy in the following form: 



E = NTu 
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Li a (e»/ T ) 
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Although Eqs. (|3.1j) and (|3.2p are exact, they are not 
very useful for practical purposes. One should consider 
them in various limiting cases, in order to obtain explicit 
expressions for fi and other thermodynamic functions. 
As will be clear, the temperature behavior of /k(T) is 
different at different values of a. If a > 2 (that is, if 
7 < 2/3, 7 < 1, and 7 < 6 at D = 1, D = 2, and 
D = 3, respectively) then one may use power expansion 
in fugacity e M / T : 
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Since /x -C T, which is valid in the vicinity of To, then, 
to the first order, e^ n l T — 1 sa fxn/T, and the summation 
is obtained in the following closed form: 



ty J_ = ii C ( a _i). 



The explicit equation for finding the chemical potential 
valid at a > 2 is 



III. THERMODYNAMIC FUNCTIONS NEAR 
BEC TRANSITION TEMPERATURE 
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The temperature behavior of the thermodynamic func- 
tions in the vicinity of the critical temperature T is found 
from the behavior of fJ,(T). Since /i(T) = at T < Tq, one 
should consider the region T = T +AT, where AT < T . 
H(T) is found from Eq. (gTFJ) : 



t; Li a (i) 



(3-1) 



where Li Q (z) is the polylogarithm function. Its defini- 
tion and some properties necessary in the context of the 
present work, are gathered in Appendix. With the help 
of Eq. (jAll) in the Appendix one can find the internal 
energy in terms of fi = /x(T). To this end one should use 
the following relations: 

E = Aj o e( jJ;_ i = AT^T(a + l)U a+1 (e^) 

a T(a)U a (>/ T ) . 
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whose solution at T > To is 
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One can see from this expression that n(T) is continuous 
at Tq, but its first derivative has a discontinuity 
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Hereafter 



[/] = lim e ^ [/(T + e) - /(T - e)} 

designates the discontinuity of the function /(T) at the 
BEC transition temperature T . Hence, the energy of the 
Bose gas is continuous, but the heat capacity at constant 
external field defined as 



is discontinuous at T = Tq, i.e., 

e a e e / T de 



The case when the inequality 
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1 < a < 2 



(3.4) 



(3.5) 



is satisfied, demands the separate treatment, because 
e /W T unc j er the sum m gq p 3[) cannot be expanded in 
powers of /x/T due to the divergence of the summation 
over n. In due turn, both the cases of the strict inequal- 
ity 1 < a < 2 and the equality a = 2 in (I3.5|) should 
also be treated separately. Depending on the number of 
space dimensions D = 1, 2, and 3, the strict inequality 
is satisfied for the power 7 in the intervals 2/3 < 7 < 2, 
7 > 2, and I7I > 6, respectively. 

I 



First, let us consider the case of the strict inequality 
1 < a < 2. At small /1, the dominant contribution to the 
sum in Eq. (I3.3[) comes from large n, so the summation 
can be approximated by the integration with the help of 
the formula due to Euler and Maclaurin Q: 

E /(« + «)« / f(x)dx - -f(a) - -/'(a).(3.6) 

n=l Ja 

To this end one should write 

00 00 00 

E = E /(* + 1) = ^) + E /(* + !). 



z=i 



and apply (|3.6p in the particular case a = 0. Upon ne- 
glecting the non-singular terms of the order of fi/T, one 
finds that the dominant contribution to the equation re- 
lating the chemical potential fi to the particle number N 
is thus represented in the form 
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where 
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is the incomplete gamma function Q, and a; = |/x|/T <C 
1. Because x is small, one can keep only the term with 
n = in the sum over n in Eq. (|3.8|) . As a result, after 
keeping the terms of the lower order in the ratio fi/T, one 
obtains the first-order expression for the nonzero value of 
/i(T) for T > T : 



a(a-l)C(a) (T-T ) 



T r (2 - a) 



(3.9) 



This expression shows that the chemical potential is con- 
tinuous at T = Tq. Indeed, (J>(T) vanishes at T < T 
with all its derivatives. Since, in the present case, 
l/(a-l) > 0, then, at T > T , the expression (T-T )^ 
vanishes at T —> Tq, too. As for the derivatives of the 
chemical potential over temperature, the evaluation of 
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the order k derivative can be expressed as follows 



d k /i 
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which shows that, provided the condition 



1 



(3.10) 



is fulfilled, where k > 2, the A:-th order derivative of the 
chemical potential as a function of temperature is discon- 
tinuous at T = Tq. The magnitude of the discontinuity 
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(3.11) 



Since, in the present case, a < 2, the case k = 1 is beyond 
the treatment. As will be clear later on, when a = 2, the 
temperature behavior of the chemical potential becomes 
non- analytical at T = Tq. 
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Let us turn to the case a = 2. Note that this equality 
takes place, if, at the space dimension D = 1, 2, 3, the 
power of the coordinate dependence of the trapping po- 
tential U(r) oc r 1 equals, respectively, 7 = 2/3, 2, 6. It is 
the situation D = 2, 7 = 2 which was realized recently in 
the experiments with the Bose - Einstein condensation 
of photons in the cavity Q. The chemical potential at 
T > T Q can be found from Eq. (pTTj) . with a = 2: 



Li 2 (e"/ T ) Li 2 (1 - M/T) 



Li 2 (l) 



Li 2 (l) 



(3.12) 



The corresponding special function is now Li 2 (x). Its 
definition and some properties necessary for the present 
treatment, are listed in the Appendix. Using Eq. (|A3[) 
with e = \fi\/T, the equation for the determination of \fjt\ 
at < T- T Q < T , reads 



T T 



n 2 (T-T ) 
3T : 



whose solution, with the logarithmic accuracy, can be 
represented in the form 



n 2 (T-T ) 



3 In 



„ 3T 

k 2 (T-T ) 



(3.13) 



One can see that, in the experimentally accessible case 
D = 2, 7 = 2 5], the temperature dependence //(T) 
is non- analytical at T = To. The differentiation over 
temperature shows that the chemical potential and its 
first derivative are continuous but the second derivative 
has an infinite discontinuity at the critical temperature. 
Hence the temperature derivative of the heat capacity is 
discontinuous. 

Let us comment on the the thermodynamics of the 
Bose gas in the limiting cases 7 = or 7 — > 00. The first 
one is reduced to the constant external potential U(r) — 
Uq =const. In this case, the thermodynamical properties 
are meaningful, if one restricts the system to the finite 
D-dimensional volume. Since the effect of the potential 
U (r) =const can be removed by the change of the energy 
reference point, the thermodynamical properties at 7 = 
are the same as in the case of free bosons. Curiously, but 
the seemingly opposite case 7 — > 00 corresponds exactly 
to the same case of the free Bose gas. The only difference 
is that the effect of the finite volume is reached naturally. 
See Eq. $n$. The behavior of the BEC fraction N (T) 
and the chemical potential /i(T) near the BEC transition 
point To can be inferred from the results obtained in the 
preceding section by making the replacement a — > D/2. 
In particular, 



iVo(T) = N 



D/2 



which coincides with the results of [4( . For the ideal Bose 
gas, BEC is possible when D > 2. If D > 4, the chemical 



potential is continuous at T = To, but its derivative over 
T has the discontinuity of the magnitude 



[A = 



2CP/2-1)' 



If D/2 = 1 + l/k, where k = 2,3,---, that is, D = 
3, 8/3, then the second, third, ■•• derivative of the 
chemical potential is discontinuous. In fact, one finds 
from (|3.11[) , that, in the three-dimensional space, the 
magnitude of the discontinuity is 



[/<" 



9C 2 (3/2) 
8^T n 



2.44 



This coincides with the textbook result. Finally, if D = 4, 
then fi(T) for T > To is given by the expression (|3.13p . 
the same as in the case of the oscillatory external poten- 
tial in the two-dimensional space. 



IV. DISCUSSION AND CONCLUSION 



The temperature behavior of the chemical potential as 
the function of the temperature is important from the 
point of view of establishing the kind of the phase tran- 
sition, if such a transition indeed takes place. There is 
a well-known Ehrenfest classification of the phase tran- 
sitions according to which the transition is of the first 
order, if the energy (in particular, the chemical poten- 
tial) is discontinuous, of the second order, if the deriva- 
tive of the energy over temperature (the heat capacity) 
is discontinuous, etc. If one treats the Bose - Einstein 
condensation as the phase transition, then the results of 
the preceding section help to establish the kind of this 
transition. One can see, that in all cases except the one 
with a = 2, the BEC transition is indeed the phase tran- 
sition of some particular kind. Specifically, when a > 2, 
the phase transition is of the second kind [see Eq. (|3.4j) ] , 
when a = 1 + l/k, where k > 2, the transition is of 
the (k + l)-th kind [sec Eq. (|3.11|) ]. But in the singular 
case of a = 2, which is realized in D = 4 for the bosons 
in zero external potential (unphysical case), and in the 
two-dimensional oscillator (physically realized in Q), the 
BEC transition is the phase transition of the third kind. 

The kind of the phase transition could be established 
in the caloric experiments, i.e. via the measurement of 
the temperature dependencies of various thermodynam- 
ical quantities. Is the crossover character (|3 . 13[) of the 
BEC transition realized or not, was not studied in the 
experiment Q with the D = 2 Bose gas of photons in the 
oscillator potential U(r) cx r 2 . The fact is that the Bose 
- Einstein condensation in the experiment [5j took place 
at the fixed room temperature. BEC was manifested as 
the macroscopic occupation of the state with the energy 
E n = hcs/D , see Eq. ([l.ip . upon reaching some critical 
pump power of the laser corresponding to N ~ 10 4 pho- 
tons in the cavity. This type of the experiment is not the 
caloric one. 
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To summarize, the power 7 of the trapping spherically 
symmetric potential, U oc r 7 , for dimension of the space 
D which admits the Bose - Einstein condensation of the 
ideal Bose gas, is found. The degree of singularity of the 
chemical potential and other thermodynamic functions 
appears to be very different for different space dimen- 
sions D and different power 7 of the coordinate depen- 
dence of the trapping potential encoded in the quantity 
a = D/2 + D/7. In particular, the case a = 2 cor- 
responding to D = 2 and 7 = 2, accessible through 
the recent experiment || with photons in the cavity, 
results in the non-analytic behavior of the temperature 
dependence of the chemical potential and other thermo- 
dynamic functions. The results obtained in the paper, 
could be possibly applied in the case of the hypothetical 
axionic BEC as the dark matter [||, in some power-like 
spherical-symmetric potentials in 3D, realized in the self- 
gravitating distributions of the above hypothetical parti- 
cles. 

The author thanks an anonymous referee for the care- 
ful refereeing and the constructive remarks and sugges- 
tions. I am grateful to Prof. Martin Weitz for correspon- 
dence which helped to reveal the error in the first draft 
of the manuscript. 
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have been informed by Prof. V. I. Yukalov about his 
results [l2] on similar subject. 

Appendix A: Polylogarithm function 



The polylogarithm function Li Q (z) is defined by the 
following integral and series representations fTpj ] : 

_ poo a-lj 00 k 

L U Z ) = ^- - — ^ = y£_ (A1) 
a{ ' T(a)J e*-z ^k<* 1 ' 

This function is defined at \z\ < 1, and a > 1. As is 
clear from this definition, Li Q (l) = C( a )- I n particular, 



another equivalent representations for Li2(z) is used in 
the paper: 

f z dt 
Li 3 (z) = - / bi(l-t)-. (A2) 
Jo t 

To prove this representation, one should use the series 
expansion of ln(l — t) [ll|. Then 

ln(l -t) _ N " t^ 1 
t k ' 

and 

r z fh °° -? k 

-[ Mi -i)f-E i,-U,M. 

J{> fe=i 

The expansion of 1^2(2:) near z = 1 used in the body of 
the paper, is 

Li 2 (l-e)« ^r+elne, (A3) 

where e -C 1. Indeed, one has the following chain of 
equations: 

Li 2 (l-e) = -^ £ ln(l-i)y- 

U 2 (l)+j\n(l-t)j = 

L12 1 + / « 

Jo 1 - u 

— + elne. 
6 

The approximate equality has the logarithmic accuracy, 
and the relation £(2) = 7r 2 /6 is taken into account. 
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